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Abstract. We present the results of realistic, 3D, hydrodynamical, simulations of surface 
convection in red giant stars with varying effective temperatures and metallicities. We use 
the convection simulations as time-dependent, hydrodynamical, model atmospheres to com- 
pute spectral line profiles for a number of ions and molecules under the assumption of local 
(~| \ thermodynamic equilibrium (LTE). We compare the results with the predictions of line for- 

Oj. mation calculations based on ID, hydrostatic, model stellar atmospheres in order to estimate 

I ' the impact of 3D models on the derivation of elemental abundances. We find large negative 

O ' 3D- ID LTE abundance corrections (typically -0.5 to - 1 dex) for weak low-excitation lines 

from molecules and neutral species in the very low metallicity cases. Finally, we discuss the 

^ ■ extent of departures from LTE in the case of neutral iron spectral line formation. 

C^ . 
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^^^ ' 1. Introduction dimensional nature of convection, the applica- 



tion of time-independent hydrostatic ID model 
Classical spectroscopic abundance analyses of atmospheres to spectroscopic analyses is there- 
late-type stars normally make use of theoret- fore a potential source of severe systematic er- 
ical ID model atmospheres constructed un- rors. Moreover, analyses based on ID models 
der the assumptions of plane-parallel geome- cannot predict strengths and profiles of spec- 
try or spherical symmetry, hydrostatic equilib- tral lines without invoking ad-hoc fudge pa- 
rium, and flux constancy. In late-type stars, the rameters (e.g. micro- or macro-turbulence). In 
convection zone reaches and affects the stellar recent years, on the other hand, it has be- 
surface layers from which the outgoing stel- come possible to perform realistic 3D hydrody- 
C^ ■ lar flux emerges; one-dimensional model stel- namical simulations of stellar surface convec- 
lar atmospheres, however, inherently rely on tion (e.g. Nordlund 1982; Nordlund & Dravins 
rather crude implementations of convective en- 11990; Stein & Nordlund 1998; Asplund et al. 
ergy transport such as the mixing length the- M999; Frevt ag et al.ii20 02; Ludwig et al.ii2002t 
cry (lB6hm-Vitenselll958l) or analogous alter- Icarl sson et al.l l2004t 1V6gler 2004). In the case 
native formulations (e.g. 'Canut o & MazzitelUl of solar-like stars, such 3D simulations have 
Il99lh . In view of the dynamic and multi- 
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been proven successful in reproducing self- 
consistently observational constraints such as 
the topology of the granulation pattern and 
the detailed shape of spectral lines (see re- 
view by Asplund 2005) . Here, we will present 
some results of our recent 3D simulations of 
convectio n at the surface of red giant stars 
dCoUet et al. 2006, 2007) focusing in particular 
on their application to spectral line formation. 



2. The hydrodynamical simulations 

W e use the 3D radiat i on-hy drodynamics code 
of IStein & Nordlundl ( Il998l) to simulate con- 
vection at the surface of red giants with varying 
effective temperatures (reff*4700 . . .5100 K), 
a surface gravity of log g - 2.2 (cgs), and dif- 
ferent metallicities ([Fe/H]= 0... - 3). The 
radiative-hydrodynamical equations are solved 
on a discrete mesh (100x100x125) for a rep- 
resentative volume of stellar surface covering 
about eleven pressure scales in depth and at 
least ten granules horizontally, and extending 
from log Tgooo^ - 4 to logT5ooo~7 in terms of 
continuous optical depth at /I = 5000 A. Open 
boundaries are employed at the top and bot- 
tom of the domain and periodic boundaries 
horizontally. The simulat ions make use of re- 
alistic equatio n-of-state JMihalas et al. 19881) 
and opaciti es JGustafsson et al.l [19751; iKurucJ 
Il992l fl993'). We adopt the solar chemical com- 



position from Anders & Grevesse (1989) with 
the abundances of all metals scaled proportion- 
ally to [Fe/H]. At each time-step, the radiative 
transfer equation is solved along the vertical 
and eight inclined rays; opacities are grouped 
into four opacity bins (Nordlund 1982) and lo- 
cal thermodynamic equilibrium (LTE) without 
scattering terms in the source function is as- 
sumed. 

Qualitatively, the atmospheric structures 
and gas flows predicted by the 3D red gi- 
ant convection simulations are similar to 
the ones previously reported for solar-type 
stars (Stein & Nordlund 1998; Asplund et al. 
ll999HAsplund & Garcia Perezll200lh . Figured 
shows the atmospheric temperature struc- 
ture of two simulation snapshots at solar 
and [Fe/H]= -3 metallicity; the stratifica- 
tion from ID, plane-parallel, hydrostatic marcs 



model atmospheres (iGustafsson et al.l Il975t 
lAsplund et al.l fl997h generated for identical 
stellar parameters, input data, and chemical 
compositions are also displayed. At metallic- 
ities near solar, the mean 3D thermal stratifi- 
cation of the upper photosphere remains close 
to the ID structure where radiative equilibrium 
is enforced; at very low metallicities on the 
other hand, the surface layers of 3D simula- 
tions tend to be substantially cooler than in ID 
models. In 3D, the temperature in the upper 
photosphere is primarily regulated by the com- 
petition between adiabatic cooling due to the 
mechanical expansion of the gas and radiative 
heating due to reabsorption by spectral lines of 
radiation emitted from deeper in. At very low 
metallicity, because of the scarcity and weak- 
ness of spectral lines, adiabatic cooling domi- 
nates, hence the balance between cooling and 
heating is reached at lower temperatures than 
in stationary ID models (.Asplund et al , 1999.). 

3. Spectral line formation 

We use the red giant simulations as 3D, hy- 
drodynamical, time-dependent, model atmo- 
spheres to calculate spectral line profiles of 
various ions and molecules under the assump- 
tion of LTE. We compute flux profiles for about 
80 wavelength points per spectral line using the 
same radiative transfer solver as in the convec- 
tion simulations. We stress that, contrary to ID 
calculations, 3D line formation does not rely 
upon free parameters such as micro-or macro- 
turbulence to account for non-thermal Doppler 
broadening and wavelength shifts due to bulk 
motions of the gas in the photosphere. 

We estimate the impact of 3D hydrody- 
namical models on spectroscopic analyses, by 
differentially comparing the curves-of-gwwth 
of spectral lines computed with 3D and ID 
MARCs models using exactly the same input 
physics and numerical schemeQ Systematic 
differences between the 3D and ID ther- 
mal structures as well as velocity gradients 
and temperature and density inhomogeneities 
present in the 3D models, can translate into sig- 



' For the ID cases, we consider two choices of 
micro-turbulence, ^ = 1.5 kms^'and ^ = 2.0 kms"'. 
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Fig. 1. Atmospheric structures of two 3D red giant simulations at solar (left) and [Fe/H]-3 metallicity 
(right). Gray shaded area: temperature distribution as a function of continuum optical depth atA = 5000 A; 
darker areas indicate values with higher probability. Dashed line: mean temperature stratification (i.e aver- 
aged temperature on surfaces of constant optical depth) of the 3D model atmospheres. Solid line: tempera- 
ture stratification of corresponding ID marcs model atmosphere. 



nificant differences in terms of line strengths 
predicted by the two kinds of models for a 
given chemical composition. This also im- 
plies that the results of spectroscopic abun- 
dance determinations will, in general, depend 
on whether ID or 3D models enter the analy- 
sis. 

4. Results 

Figure |2] shows the differential 3D- ID LTE 
Fe abundances derived from Fe i lines at /I = 
5000 A with varying strengths and excitation 
potentials for two red giants at solar and very 
low metallicity. In LTE, at a given Fe abun- 
dance, weak Fe i lines tend to be stronger 
with 3D models than they do in ID, implying 
that negative 3D-1D LTE abundance correc- 
tions are expected. Such coiTections are more 
pronounced at low metallicities and for low- 
excitation lines: the differential 3D-1D LTE 
Fe abundances derived from weak Fe i lines for 
the model at solar metallicity are ^ -0.1 dex, 
while they can be as large as -0.8 dex at 
[Fe/H] = -3. As line strength increases, differ- 
ential 3D-1D LTE Fe abundances remain large 
and negative in the very metal-poor case, while 
they grow positive at solar metallicity. We note, 
however, that the trend in the corrections for 
strong lines at solar metallicity might well be 
an indication that the present simulations do 



not resolve the full spectrum of convective ve- 
locities and therefore underestimate the non- 
thermal Doppler broadening. Preliminary test 
simulations of surface convection in a red giant 
at [Fe/H]= and at a numerical resolution of 
200x200x250 show indeed a more prominent 
high-end tail in the distribution of convective 
velocities than at lower resolution. 

Qualitatively, we can interpret the be- 
haviour of differential 3D-1D Fe abundance 
corrections by looking at the variations of the 
number density of neutral Fe with depth in ID 
and 3D model photospheres (Fig. [3]l. At so- 
lar metallicity, temperature and density vary 
with optical depth in a fairly similar way in 
both the ID and mean 3D stratifications and 
so does the fraction of neutral to total Fe. 
Consequently, the 3D-1D abundance correc- 
tions at solar metallicity have to be ascribed 
to temperature and density inhomogeneities 
and velocity gradients in the 3D models rather 
than to differences between the ID and mean 
3D structures. The situation is radically differ- 
ent at very low metallicities. In the ID marcs 
stratification, iron is nearly completely ion- 
ized throughout the stellar atmosphere; in the 
upper photospheric layers of the metal-poor 
3D model atmosphere instead, we observe a 
bimodality in the distribution of the fraction 
of neutral Fe versus optical depth. While in 
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Fig. 2. 3D- ID LTE Fe abundance corrections derived for two red giants at solar (left panel) and at [Fe/H]= 
-3 metallicity (right panel) from Fe i lines a.lA = 5000 A as a function of equivalent width W^. Corrections 
are shown as a function of lower level excitation potentials of the lines and for two different choices of 
micro- turbulence in the ID calculations. 
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Fig. 3. Ratio of neutral to total Fe number densities («pe,/«Fe) as a function of optical depth aiA = 5000 A 
in the atmosphere of two red giants, at solar (left panel) and at [Fe/H]= -3 metallicity (right panel). The 
gray shaded area represents the distribution of np^j/np^ values with optical depth in the 3D model atmo- 
sphere; darker areas indicate values with higher probability. Over-plotted are the curves for the mean 3D 
stratifications (dashed line) and the corresponding ID marcs model atmospheres (solid line). 



the hottest parts of the upper 3D photosphere 
iron is nearly completely ionized, in the at- 
mospheric regions where the temperature falls 
below a critical value (~3000 K) iron is pre- 
dominantly neutral. In other words, due to the 
systematically cooler temperatures in the up- 
per photosphere of 3D metal-poor models, Fe i 
lines tend to be much stronger in 3D than in 
ID, leading to large negative 3D-1D LTE cor- 
rections to the Fe abundance. From a qualita- 
tive point of view, neutral lines of other ele- 
ments or molecular features behave similarly 
to Fe I lines and are also characterized by 



large 3D-1D LTE abundance coiTections at 
very low metallicity. The actual value of the 
corrections depend on the details of the ioniza- 
tion and molecular equilibria. At this stage, it 
is important to caution that, Fe i lines are most 
likely affected by departures from LTE. The 
main non-LTE mechanism for Fe i in late-type 
stellar atmospheres is efficient over-ionization 
driven by the UV radiation field. This causes 
an underpopulation of Fe i levels with respect 
to LTE, and leads to weaker Fe i lines, and, in 
turn, to the prediction of higher Fe abundances 
in non-LTE. ICollet et al] (l2006h have estimated 
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Table 1. Fe abundance derived for the red gi- 
ant HE 0107-5240 using a ID marcs model 
and the mean stratification from a 3D model. 
The non-LT E calculati ons assume the model 
Fe atom from Collet et al. (2005) . with efficient 
inelastic H+Fe collisions and thermaUzation of 
the uppermostFei levels. 



ther details see lCoUet eFani2007h . Finally, we 
have also examined possible departures of Fe i 
line formation from LTE; such non-LTE cor- 
rections are opposite to and, according to pre- 
liminary ID test calculations of ours, of the 
same order of magnitude as the ones due to 
granulation. 



Model [Fe/H]LTE [FeAi]„o„-LTE 



ID 

mean 3D 



-5.40 
-5.60 



-4.65 

-4.75 



the non-LTE effects on Fei lines for the ex- 
tremely iron-poor giant HE 0107-5240 using 
both a ID MARCS model atmosphere and a mean 
atmospheric stratification from 3D simulations 
for the analysis. The estimated non-LTE effects 
(TablelTJ are considerable and opposite to the 
3D-1D LTE corrections, implying that a com- 
bined treatment of 3D and non-LTE eff'ects is 
crucial for accurate Fe abundance determina- 
tions. 



5. Conclusions 

We have presented here some illustrative re- 
sults of the application of 3D red giant sur- 
face convection simulations to spectral line for- 
mation in LTE. The differences between the 
temperature stratifications of 3D simulations 
and corresponding ID model atmospheres as 
well as the 3D temperature and density inho- 
mogeneities and velocity gradients can signif- 
icantly affect the predicted line strengths and, 
consequently, the value of elemental abun- 
dances inferred from spectral lines. At very 
low metallicities, where the deviations of the 
mean 3D thermal structure from the classical 
ID stratification are largest, the 3D- ID LTE 
abundance coiTections are negative and con- 
siderable for lines of neutral species (about 
-0.8 dex for weak low-excitation Fei lines). 
Corrections to CNO abundances derived from 
CH, NH, and OH weak low-excitation are also 
found to be typically in the range -0.5 dex to 
-1.0 dex in very metal-poor giants (for fur- 
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